Thermal Relaxation Time in Chemically Non-equilibrated 

Quark- Gluon Plasma 

Xiao-Fei Zhang Wei-Qin Chao'^'* 
"CCAST (World Laboratory), P. O. Box 8730, Beijing 100080, People's Republic of China 
^Institute of High Energy Physics, Academia Sinica, Beijing 100039, People's Republic of China 



Abstract 

The definition of thermal relaxation time is extended to chemi- 
cally non-equilibrated quark-gluon plasma and the chemical non- 
equilibrated thermal relaxation times for partons are calculated 
using the non-equilibrium Debyc mass as the infrared regulator. 
The dependence of the thermal relaxation time on the fugacity 
is given and the influence of the chemical non-equilibration is 
discussed. We find that there are threshold fugacities A* and 
a; for gluons and quarks, for Xg < A;(A, < X*), t^'^^/t^^ 
(^^NEQ i^EQ \^ decreases strongly with increasing fugacity, while for 
■^*g{\ > Ki)^ ratios are almost 1. It is shown that there 
is also the two-stage equilibration in a chemically non-equilibrated 
plasma. We also discussed the effect of using the non-equilibrium 
Deby mass as the infrared cutoff. 



PACS number(s): 12.38.Mh, 24.85.+p, 52.25.Dg 



I. INTRODUCTION 

One of the main objectives of the future experimental programs at the BNL Rel- 
ativistic Heavy Ion Colhder (RHIC) and the CERN Large Hadron Collider (LHC) is 
the production of Quark-Gluon Plasma (QGP) through heavy ion collisions [1]. The 
parton gas produced in heavy ion collisions is not immediately in thermal and chem- 
ical equilibrium. The parton gas relaxes to the equilibrium state through secondary 
interactions. Because the preequilibrium phase may influence the QGP signals, it 
is essential to discuss nonequilibrium properties of QGP and to investigate whether 
QGP produced in RHIC and LHC could reach equilibrium. Quantum transport 
theory[2] in principle can be used to describe the properties of nonequilibrium QGP 
and there are some discussions based on it [3]. However, its applicability to heavy 
ion collisions is still far from being realistic. The Boltzmann equation in the re- 
laxation time approximation is often used to discuss the transport coefficients and 
the evolution of the QGP [4-5]. In this approximation method the collision terms 
are determined by the corresponding thermal relaxation time which measures the 
time scale of the nonequilibrium system approaching equilibrium and can be esti- 
mated from the parton interaction rates [6-8]. The equilibration of QGP considering 
of parton expansion effect is discussed recently [5] However, the thermal relaxation 
time has been considered so far mostly for a chemically equilibrated plasma. As 
discussed by many models[9-ll], in the early stage of ultra-relativistic heavy ion 
collisions, the parton gas is dominated by gluons and far from chemical equilibrium. 
First, the system is cvoluting to thermal equilibrium mainly through clastic scat- 
tering, which leads to the local isotropy in momentum distribution. The system 
evolute towards chemical equilibrium later and sometimes the plasma could never 
reach chemical equilibrium during its life time. In this paper we calculate the ther- 
mal relaxation time in a chemically non-equilibrated plasma and discuss how the 
chemical nonequilibration influences the parton thermal relaxation. 

The thermal relaxation time in QGP diverges strongly at small momentum trans- 
fer in the naive perturbation theory. This difficulty could be overcome by the 
Bratten-Yuan method[12] which proposed that the soft contribution k < k* and 
the hard part k > k* could be calculated separately by introducing a separation 
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scale k* for the momentum transfer. The soft part is treated using the resumed 
propagators and vertices proposed by Braaten and Pisarski[13], whereas the bare 
Green functions are sufficient for the hard one. Assuming gT < k* < T, the final 
result is independent of the separation k* . There is another much simpler and widely 
applied approximation method based on the naive perturbation theory, using bare 
propagators with the Debye mass as infrared regulator [14]. It was demonstrated 
that this method also worked well for calculating the thermal relaxation time[6]. 
We will use this approximation approach, and the nonequilibrium Debye mass as 
the infrared regulator in this paper. 

This paper is organized as follows, in Section II we first define the interaction 
rates in a chemically non-equilibrated plasma. Then, in Section III, the thermal 
relaxation time for quarks and gluons are calculated using the definition of section 
II. Finally, in section IV, we give some numerical results and discussions. 



II. DEFINITION OF THE INTERACTION RATES IN A 
CHEMICALLY NON-EQUILIBRATED GAS 



In this section we will extend Wcldon's discussion [17] and define the interaction 
rates for a chemically non-equilibrated gas. There are two kinds of interaction rates. 
One is the ordinary interaction rate defined as the damping rate for partons, which 
is usually called interaction rate. The other is the so called transport interaction 
rate obtained from the ordinary one by introducing a transport weight containing 
the scattering angle in the center of mass system. Because large angle scattering 
is the most efficient mechanism for the dissipative momentum transfer in the case 
of a plasma with long range interaction, hence the mean free path, as well as the 
thermal relaxation time, may rather be defined by the inverse of the transport 
interaction rate [15]. Let us first define the ordinary interaction rate for a chemically 
non-equilibrated plasma from the Bolzmann equation[16]. 




[hh{l ± /3)(1 ± U) - (1 ± ± h)hh\ 
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(27r)^5^(Pi + P2-P3- P4)|Mi2^34r, (1) 

where Pi — ipi,Pi), = P^/Pi- 5f = d/d{xi)fj,. fi is the distribution function f{pi) 
for partons. ± are used for bosons (gluons ) and fermions(quarks and anti-quarks). 
I -^12-^34 P is the squared matrix element for corresponding parton scattering process 
and is summed over final states and averaged over initial states. The spin and color 
factor 1/ is 16 for gluons and 6 Nf with Nf flavors of quarks or anti-quarks. 

To obtain the definition of the interaction rates [17], we assume that /i is away 
from equilibrium and fi,i = 2,3,4, arc in thermal equilibrium but not in chemical 
equilibrium, which can be expressed as the following form, 

// = A.^^^, i = 2,3,4, (2) 

where Aj is the fugacity of the corresponding partons. It measures how far the system 
is from the chemical equilibrium. Using the relation for the distribution function 
Eq.(2), 

l±/f = K'e^^-n, (3) 
we obtain from the Boltzmann equation, 

^^i-ai/i(pi) = -r(pi)(/i-/r), (4) 

{27Ty6\P^ + P2-P3- P4)|Mi2^34r. (5) 

The above equation is the definition of the interaction rates in a chemically non- 
equilibrated plasma. By putting Aj = 1, Eq.(5) is just the interaction rates for a 
chemically equilibrated plasma[17]. 

We will assume = /|,/| = fl in the following discussion. This simplification 
holds as long as k = \pi — p^l = \p2 — Pi\ is not too large or —t is not of the order 
of s, with s,t,u being the usual Mandelstam variables [8, 16]. Now the transport 
interaction rates and the thermal relaxation time can be defined from Eq.(5) after 
introducing a transport weight sin9^/2, 



r d^P2 


f d^Ps 


f d^p, /|(1±/|)(1±/|) 


1 (27r)32p2^ 


(27r)32p3 J 


(27r)32p4 (l±/f) 
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T 2pJ (27r)32p2i {2tt)^2psJ {2tt)^2p/^^ 

{2ti)H\P, + P2-Ps- P4)|Mi2-.34r^. (6) 



III. RELAXATION TIME FOR QUARKS AND GLUONS 

In this section we use Eq. (6) to calculate the thermal relaxation time for quarks 
and gluons in chemically non-equilibrated plasma. As thermahzation is achieved 
mainly through momentum changes in elastic scattering to the leading order. Here 
we only consider the contribution of the elastic scattering. In the lowest order the 
relevant matrix elements are [18] 

K-.-r^^'^. (s) 

M = . (9) 

Using the definition of the differential cross section [19] 

da 



dt 167rs2 

for gluons, we obtain 



MY (10) 



+12JV/ / /'(*)[! - /,'(*)] / *(f (11) 

Where fg and fq are the chemical non-equilibrated distribution function of gluons 
and quarks respectively. 

The interaction rates diverges at small k. As has been discussed, for a chemically 
equilibrated plasma, it can be regulated by the Debye mass which is the gluon 
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self energy at zero momentum in the high temperature hmit[20]. For a chemically 
non-equilibrated plasma, we use the non-equilibrium Debye mass as the infrared 
regulator, 

ml^EQ = MK + -^KhsT'. (12) 

In the following calculation, we assume fi can be approximated by its factorized 
form[10]. 

Using Eqs. (11-13) and after some direct calculations, we obtain the thermal relax- 
ation time for gluons 

1 79T^ -K^ 

-al\2\m - 2AV(3) + Aj-] \ln{sjcf^<^) - 19/15] 



,+^^«?[iM(3) - \\\m)^y?^\\ln{sj^^^^) - 1.26], (14) 

where ^(3) = 1.202, the infrared cutoff q^^^ — rnjjj^EQ. We assume that SggUndSgg 
in the above equation can be replaced by its thermal average value, < Sgg >— 2 < 
Pg >< Pg >,< Sgg >— 2 < Pg X Pq > whcrc the expression for < > is 

<Pi>= — —k ■ (15) 

/(I^/(a) 

It can be seen from Eq(15) that < s > is approximately the same for the chemically 
non-equilibrated and equilibrated gluon gases i.e., < Sgg >= 14.59T^,< Sgg >= 
16.96T2. 

For the thermal relaxation time of quarks we consider the processes qq — > qq, 
qq — > qq and qg qg. The thermal relaxation time for quarks we obtained is 

1 32T2 tt^ 

al[2XgaS) - 2Aje(3) + A5-][/n(s,,/g^^«) - 19/15] 



Tn T^^gq ^ ^3 



'+^'irf^"'[^M(3) - ^AV(3) + X'-][ln{Sgg/q''^Q) - 1.26] (16) 



37rsqg ''2 ' 2 ' ^ 6 
Where Sqq — 19.72T^ in the above equation. 
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IV. NUMERICAL RESULTS AND DISCUSSIONS 

The thermal relaxation time depends on the fugacity of gluons and quarks. To 
completely determine it we must know the initial conditions of heavy ion colli- 
sions, which only can be obtained from models now. Focusing on the effect of the 
chemical non-cquilibration on the thermal relaxation time in this paper, we give 
a general discussion here. Fig. 1-2 show the ratio of the thermal relaxation time 
in chemically nonequilibrium QGP to the one in chemical equilibrium, Tf^Eql'^EQ 
and t%eq/'^eqi^-^- fugacity A^, Ag respectively. It is found that the results depend 
weakly on the ratio between the gluon fugacity and quark fugacity. We give the 
results of two cases, A^ = A^ and \ — \\g. j^^Q jj^Q decreases as the fugacity 
increases in both cases. There are corresponding threshold fugacities A* and A*, For 
Ag < A* (Aq < A*), t^^^/t^^ (r^^^/rq^) depends strongly on the fugacity, while 
for Xg > A* {Xg > A*), the ratio is almost 1. These threshold fugacities depend 
on the QCD coupling constant. From Fig. 1-2, one can also see that the bigger 
the coupling constant, the smaller the threshold fugacity. The results of Fig. 1-2 
can be explained as follows: When the system is approaching chemical equilibrium, 
more and more partons are produced and the fugacity increases continuously. The 
parton density of chemically non-equilibrated plasma is smaller than the one in the 
equilibrium state at the same temperature. As a result, the relaxation time in the 
chemically non-equilibrated plasma is longer than the one in chemical equilibrium. 
We also check the influences of using non-equilibrium Debye mass as the infrared 
cutoff and the result is shown Fig. 3. Using the nonequilibrium Debye mass as 
the infrared regulator restrains the increase of r^^Q/r^Q with decreasing fugacity 
strongly. 

Based on parton interaction rates calculated for the chemically equilibrated 
plasma, Shuryak argued that equilibration of the plasma proceeds via two stages 
in "hot gluon scenario" [21], where gluons equilibrate much faster than quarks. To 
investigate if the situation is changed for a chemically non-equilibrated plasma we 
also calculate the ratio Tg/Tg in a chemically non-equilibrated plasma. The results 
are shown in Fig.4. We find that the ratios are around 2 and do not change much 
with the fugacity. although it changes a little more for longer coupling constants. 
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Prom our calculation we can conclude that there is also the two-stage equilibration 
in a chemically non-equilibrated plasma, i.e., gluons also equilibrated much faster 
than quarks. 

As a final remark, we should point out that only the contribution of the elas- 
tic scattering to the thermal relaxation is taken into account in this work. The 
contribution of the inelastic scattering should also be considered in the future study. 
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Figure Capure: 

Fig.l. The ratio of the thermal relaxation time of gluons in chemically nonequi- 
hbrium plasma to the one in chemical equihbrium, Tg^^ /Tg^, v.s. fugacity A^. 

Fig. 2. The ratio of the thermal relaxation time of quarks in chemically nonequi- 
librium plasma to the one in chemical equilibrium, t^^^ /t^^ , v.s. fugacity A^. 

Fig. 3. The thermal relaxation time obtained using non-equilibrium Debye mass 
as the infrared regulator comparing with the one using equilibrium Debye mass. 

Fig. 4. The ratio of the thermal relaxation time of quarks to the thermal relax- 
ation time of gluons in chemically nonequilibrium plasma, rj^^'^/rj^^'^, v.s. fugacity 
A,. 
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Fig. 4 



